arXiv: 1507.00959vl [hep-ph] 3 Jul 2015 


Asymmetric Dark Matter Stars 


Chris Kouvaris^’* and Niklas Gr0nlund Nielsen^’ ^ 

-Origins & Danish Institute for Advanced Study DIAS, 

University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark 

We study the possibility of asymmetric dark matter with self-interactions forming compact stable 
objects. We solve the Tolman-Oppenheimer-Volkoff equation and find the mass-radius relation of 
such “dark stars”, their density profile and their Chandrasekhar mass limit. We consider fermionic 
asymmetric dark matter with Yukawa-type self-interactions appropriate for solving the well known 
problems of the collisionless dark matter paradigm. We find that in several cases the relativistic 
effects are significant. 

Preprint: CP3-Origins-2015-027 DNRF90, DIAS-2015-27 


I. INTRODUCTION 

Lately, an emerging amount of issues indicates that the 
Collisionless Cold Dark Matter (CCDM) paradigm is at 
odds with astrophysical observations. The first and most 
well known issue is related to the fact that dwarf galaxies 
have a flat density core [1, 2]. Dwarf galaxies are domi¬ 
nated by dark matter (DM) and the flatness of the density 
profile in the core of the galaxy is in contradiction with 
the cuspy profiles predicted by numerical simulations of 
CCDM [3]. Numerical simulations of CCDM predict also a 
larger number of satellite galaxies in the Milky Way than 
what has been observed so far [4-6]. Although our galaxy 
might simply be a statistical fluctuation [7-9] there might 
be dim galaxies yet to be observed, it is possible that DM is 
simply not collisionless. Furthermore, another related issue 
is the so-called “too big to fail” problem [10], i.e. CCDM 
numerical simulations predict massive dwarf galaxies that 
are too big to not have visible stars and thus to not be 
observed. Although some of the aforementioned issues can 
be resolved upon assuming the existence of baryonic-DM 
interactions [11-14], another probably more attractive and 
natural possibility is the existence of DM self-interactions. 
Clearly such interactions would flatten out cuspy dwarf 
galaxy cores and they could possibly also resolve the satel¬ 
lite galaxies issues [15-18]. 

DM self-interactions have already been proposed and 
studied in different contexts [16-35]. DM numerical sim¬ 
ulations including self-interactions favor a DM-DM cross 
section per DM mass between 0.1 — 10 cm^/g. Within 
this range, DM self-interactions can solve the cusp vs core 
problem of dwarf galaxies as well as the “too big to fail” 
one. However, DM self-interactions cannot be arbitrarily 
strong. There are several constraints imposed on them. 
First of all one should make sure that DM-DM interactions 
are not sufficiently strong to destroy the ellipticity of spiral 
galaxies [36, 37] or dissociate the subclusters of the bullet 
cluster [38]. In addition, fermionic asymmetric DM with 
attractive Yukawa-type self-interactions can lead in some 
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cases to formation of destructive black holes in the interior 
of old neutron stars, thus imposing extra constraints [25]. 
Further constraints are imposed in the case where the me¬ 
diator of the DM-DM force couples to the Standard Model. 
The mediator (p can simply couple to the Standard Model 
via e.g. a Higgs portal [39-46]. In such a case one should 
make sure that (p decays before the start of the Big Bang 
Nucleosynthesis. The fact that a minimum strength be¬ 
tween baryons and (p is required for not spoiling the BBN 
predictions can lead to significant rates of DM collisions in 
underground detectors that can exclude such models [29]. 
However constraints like these can be evaded if e.g. (p cou¬ 
ples also to sterile or active neutrinos [47] . One should em¬ 
phasize that the above constraints are model dependent and 
therefore although there is no clear universal region where 
DM self-interactions are allowed, the region 0.1 — 1 cm^/g 
which accommodate the resolution of the dwarf galaxies 
problems is roughly speaking constraint free. 

If DM experiences self-interacting forces, it is possible 
to imagine that star-like compact objects can be formed. 
Whatever the mechanisms of forming such objects, one 
should make sure that they do not violate the limits im¬ 
posed by the MACHO [48] and EROS [49] experiments. 
Based on the microlensing technique, these experiments 
claimed that less than 20% of DM can be in the form of 
compact objects between the mass range 10 “^Mq < M < 
IQMq, where Mq is the solar mass. The possibility of 
stars made of DM has been studied before in the context 
of annihilating DM forming dark stars in the early uni¬ 
verse [50-52]. It has been also studied in the context of 
hybrid compact stars made of baryonic and DM [53, 54]. 
In the latter case, neutron stars and white dwarfs include a 
signihcant amount of DM in their interior modifying thus 
the equation of state of the star. Furthermore the possi¬ 
bility of black hole formation from strongly self-interacting 
components of DM was studied recently in [55]. 

In this paper we examine the possibility that asymmetric 
DM with self-interactions appropriate for solving the core 
vs cusp problem, the “satellite problem” and the “too big to 
fail problem” forms star-like compact objects. Asymmetric 
DM [56-75] has become an attractive alternative to ther¬ 
mally produced DM not only because it can relate theories 
beyond the Standard Model to DM, but because it can also 
provide a link between baryogenesis and dark-genesis. For 
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recent reviews on asymmetric DM see [76, 77] . We are going 
to assume that the self-interactions are Yukawa-type and 
can be either attractive (mediated by a scalar (/>) or repul¬ 
sive (mediated by a vector boson (p^). Upon these assump¬ 
tions, we study the stability of dark stars formed by asym¬ 
metric fermionic DM. We solve the Tolman-Oppenheimer- 
Volkoff equation and study the hydrostatic equilibrium of 
these compact objects. We hnd their density profile, the 
mass vs radius relation as well as the Chandrasekhar mass, 
i.e. the maximum mass where these objects are stable. 

The paper is organised as follows: In section II we present 
the equation of state for DM with self-interactions. In sec¬ 
tion III we show the relevant parameter phase space of 
self-interactions that solve the problems of CCDM we have 
mentioned earlier. In section IV we present the equations 
for the hydrostatic equilibrium of the dark stars. In sec¬ 
tion V we use a simplified Newtonian analysis to get a first 
understanding of the problem, while we present the full rel¬ 
ativistic results in section VI. We conclude in section VII. 

Throughout this paper we use natural units h = c = 
/cb = I, and define the Planck mass as Mp = 


II. EQUATION OF STATE 


Pkin 
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where the functions ^ and tp are defined as 


( 5 ) 


i{x) = ^x\/ 1 -I- (1 -I- 2x^) - In X -I- \/T+~x 2 j I, 


( 6 ) 


ip{x) = |xVT+^ (2x^/3 — 1) -I- In X -I- \/l -l-x^ j |. 


_ ( 7 ) 

X = pp/m^ is a measure of how relativistic the particles 
are, and gs = 2s -I- 1 is the spin multiplicity. 

The Yukawa potential between two particles is 


Vij=±a -, ( 8 ) 

Tij 

where /r is the mass of the mediator, is the separation 
between the particles and a = is the coupling to 

the dark mediator. In order to find the Yukawa potential 
energy of the entire system we have in principle to sum 
over all pairs of DM particles, which we approximate as an 
integration over volume elements. 


As we mentioned we assume that DM is of asymmet¬ 
ric type and fermionic. We also assume that DM self¬ 
interactions are mediated by (p via a Yukawa coupling of 
the form gpxx (i^i case of attractive interactions) where x 
is the DM particle and g is the Yukawa coupling constant 
or gipfiXl^X (ill case of repulsive interactions). We wish 
to obtain the equation of state of DM under the aforemen¬ 
tioned assumptions. The energy density of DM particles p 
consists of two components 

P = Pkin + PY, (1) 


where pkin, and py are the kinetic energy density, and 
Yukawa potential energy respectively. The pressure of the 
system is [78] 


P = n 


2 


dn 



( 2 ) 


where n is the number density of DM particles. The pres¬ 
sure related to the kinetic energy P^in in the non-relativistic 
(p <C m^) or relativistic (p ^ m^) limits takes the sim¬ 
ple form of a polytrope Pkin = KPkim where T equals 5/3 
or 4/3 for the non-relativistic and relativistic case respec¬ 
tively. We choose to work with the full relativistic disper¬ 
sion Pkin = y^p^ -I- . If we assume that the temperature 

of DM particles is much smaller than their Fermi energy, 
i.e. we effectively take the limit T = 0, the number density, 
kinetic energy density and the corresponding pressure are 
given by 


n = 


/ dTTp^dp 

(27r)3 Jo 


Btt^ 


( 3 ) 


Y = 


(EC, 



- dVidVj. 


(9) 


In the case where the radius of the star satisfies R ^ 1/p 
(i.e. the potential is short range), it is a reasonable approx¬ 
imation to integrate the volume up to infinity (instead of 
the volume of the star). This leads to the following Yukawa 
energy density 


Py = ± 


2T:anJ 


p- 


I Q^g« 
187r3 


( 10 ) 


This estimate for the Yukawa energy density gives us the 
final expressions for P and p: 


P = ^m'^'tp{x) ± 

gs 4 p/ \ I 

p= 


187r3 

ags 

187r3 


X®, 

(11) 

x^ 

(12) 


Since neither equation can be inverted analytically, we must 
work with two equations of state, and have an implicit re¬ 
lation between P and p. DM particles in the attractive 
scenario (corresponding to the minus sign in the Yukawa 
contribution) cannot become arbitrarily relativistic, since 
the pressure and density must be positive. The positive¬ 
ness of pressure and density give an upper bound on x. 


III. PARAMETER SPACE OF 
SELF-INTERACTIONS 


Pkin = 


gs 

(27r)3 



E{p) Anp^dp 


2 


m^C(x), 


( 4 ) 


As we mentioned earlier, DM self-interactions within the 
range o/m^ = 0.1 — 10 cm^/g can solve the problematic 
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issues of CCDM, while astrophysical constraints limit these 
interactions between aIm^ = 0.1 — 1 cm^/g. Following [27], 
in order to determine the parameter space of DM and me¬ 
diator masses (for a given coupling) that lies in the afore¬ 
mentioned range, we introduce the transfer cross section 
tj-r = /dD(l — cos6)da/dn. We use a typical value of 
uo = 10 km/s for the average velocity of DM in a dwarf 
galaxy, and we estimate the velocity averaged cross section 
as 


cr = 


e-(«/’'o)V2 

To — 2717 ^ 

(27rUg)'^/2 


(13) 


The transfer cross section for attractive Yukawa interac¬ 
tions in the classical limit rriy^v/^ 1 is [37, 80] 

' 47r 


aT = < 


Stt 

y? 

TT 


/32iog(i + r') 


if ^<10“ 


/3V (1-b 1.5/3-1-®®) 


if 10-1 < /3 < 10® 


(log /3 -b 1 - 5 log 1 /3) if /3 > 10®. 

(14) 

The corresponding cross section for repulsive Yukawa in¬ 
teractions interactions reads 


r^/32log(l + /3-2), if P<1 
I ^ (log 2/3 - log log 2^)^ , if /3>1 


(15) 


where /3 = 2apL/(rn^v^). For the opposite limit m^v/fx < 1, 
we follow [27] and approximate the Yukawa by a Hulthen 
potential with proper parameter choices. In the region 
m^y![I ^ 1 we interpolate between the two regimes. In 
Fig. 1 we show the allowed parameter space for DM and me¬ 
diator mass (in GeV) in the case of repulsive and attractive 
interactions for three distinct values a = IQ-^, IQ-®, IQ-^. 
Deep (light) blue is the region that (j/m^ = 0.1 — 1 cm^/g 
(1 — 10 cm^/g) in dwarf galaxies, solving the aforemen¬ 
tioned issues of CCDM. The red solid (dashed) line shows 
the curve where cr = 0.1 cm^/g (1 cm^/g) in the Milky 
Way. The phase space to the left of the red curve is 
excluded because the cross section is sufficiently large to 
smooth out the ellipticity of Milky Way to a degree in¬ 
consistent with observations. There is a bit of ambiguity 
regarding the value of the maximum a consistent with 
observations but it should be between 0.1 — 1 cm^/g [27]. 


where P, M, and p are the pressure, the mass and the 
density of the star at radius r. Additionally, the continuity 
of mass gives 


dM 

dr 


Anr^p. 


(17) 


The above two equations along with the equations of state 
(11) and (12) form a complete set of differential equations 
that can be solved numerically providing pressure, density 
and mass as a function of r. As mentioned earlier, in the 
absence of self-interactions, the equation of state takes a 
simple polytropic form P = Kp^ in both the relativis¬ 
tic and non-relativistic limits. This polytropic equation of 
state together with Eqs. (16) and (17) reduce to the well 
known Lane-Emden equation with index n = 3/2, 3 at the 
non-relativistic and relativistic limit respectively. However, 
as it has been seen in the case of neutron stars, the New¬ 
tonian approximation is not accurate enough and general 
relativity must be taken into account. We implement this 
by solving the Tolman-Oppenheimer-Volkoff equation to¬ 
gether with Eqs. (17), (11) and (12). 

For completeness, let us briefly review how the Tolman- 
Oppenheimer-Volkoff equation is obtained. We are seeking 
a solution of the Einstein field equation — g^,jR/2 = 
SttGT^i,, in the presence of matter with an energy momen¬ 
tum tensor of an ideal liquid of the form 


T;=diag[p,-P,-P,-P]. (18) 

If one plugs the following spherically symmetric metric into 
Einstein’s field equations 

ds^ = - r^dn^, (19) 


one finds the following set of equations that must be satis¬ 
fied 


SttGP = e-^(- + \') 

SnGp = e~^ - 

\r J 

dP _ {P + py 

dr 2 


1 

1 

/2 ’ 


( 20 ) 

( 21 ) 

( 22 ) 


Requiring that the metric reduces to the empty space 
Schwarzschild solution at the boundary of the star, pro¬ 
vides a solution for A(r) and iz(r) and Eq. (22) takes the 
final form 


IV. STELLAR HYDROSTATIC EQUILIBRIUM 


An asymmetric dark star resembles in several aspects a 
neutron star. Both types of stars produce no energy by 
fusing nuclei in their cores. Therefore there is no radiation 
pressure present. The structure of the star is determined 
by the equilibrium between the Fermi pressure of the con¬ 
stituent elements and gravity. In Newtonian dynamics the 
above condition takes the simple form 


dP _ GMp 
dr r^ ’ 


(16) 


(W 

dr 


GMp 




4:7Tr^P 


M 


[1-^] 


(23) 


This is the relativistic version of Eq. (16). We have solved 
the coupled system of Eqs. (23), (17), (11) and (12) and we 
have obtained the structure profile of asymmetric dark stars 
both for attractive and repulsive self-interactions. Rela¬ 
tivistic effects can be quite significant and therefore it is 
compulsory to use the Tolman-Oppenheimer-Volkoff equa¬ 
tion instead of the simpler Newtonian version. 
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Logic Log^o iJ 

(e) a = 10~'* attractive (f) a — repulsive 

FIG. 1: The parameter phase space for DM and mediator masses (measured in GeV) that can solve the problems of 
CCDM. Left panels correspond to attractive interactions while the right panels to repulsive ones. We show three 
different values of the Yukawa strength a. See text for details. 
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IV. 1. Hydrostatic Stability 

Solving the Tolman-Oppenheimer-Volkoff equation 
yields an equilibrium solution, which may be stable or 
unstable. We will briefly review the conditions that must 
be satisfied in order for a star to pass from stability to 
instability. We assume a constant chemical composition 
and constant entropy per DM particle. 

The total mass-energy of the star is the integrated energy 
density 


.R 

M{R) = / Airr'^pdr, (24) 

Jo 

where R is defined by P(i?) = 0. The number of constituent 
DM particles in the star is [79] 


nR 

N{R) = / 47rr^ 
Jo 



2GM{r) 

r 


- 1/2 

ndr, 


(25) 


where n is the DM number density as a function of r. A 
star of constant chemical composition and entropy per par¬ 
ticle can only pass from stability to instability with respect 
to some particular radial normal mode, at a value of the 
central density pc for which we have [79] 


^ = 0 ^=0 
dpc ’ dp, ■ 


(26) 


These conditions are satisfied if we alternatively choose to 
satisfy simultaneously the first equation above and 


d /M\ 
Wc 


= 0 . 


(27) 


Fig. 2 shows M/N (divided by m“^) and M (in units of 
the corresponding upper mass Mch) as a function of p, 
in the repulsive case with = 100 GeV, p = 10 MeV 
and a = 10“^. One can see that the points where M/N 
minimizes and M maximizes coincide. This is a generic 
feature for all the star profiles we present in Figs. 3 and 
4, i.e. the transition from stability to instability (collapse) 
takes place at the point where the mass maximizes. 

In the absence of Yukawa interactions, the critical rela¬ 
tivity parameter at transition from stability to instability 
is x, ~ 0.8 [78], which is independent of m^. Yukawa in¬ 
teractions induce dependence of x, on p and a. 


V. ANALYTIC NEWTONIAN APPROXIMATION 

We would like to estimate the upper mass limit (Chan¬ 
drasekhar mass) for asymmetric dark stars using first sim¬ 
ple Newtonian arguments. The reader interested in the full 
relativistic results may skip this section and move to the 
next one. In order to understand the different regimes the 
dark star passes through before a collapse occurs we de¬ 
rive analytical solutions using a Newtonian approximation 
for gravity and assume a constant density. Comparing this 



central density. In this example the interactions are 
repulsive, = 100 GeV, p = 10 MeV and a = 10“^. 
The central density at which the energy per particle is 
minimized coincides with the density at the 
Chandrasekhar mass in Fig. 3 (a) (marked by a blue 
circle). For this reason the equilibrium conhgurations at 
higher densities are unstable. As expected, at low 
densities M/N is simply m^. 


approach to the full relativistic treatment of the Tolman- 
Oppenheimer-Volkoff equation, will provide an idea of how 
important general relativity effects are in different DM sce¬ 
narios. 

In this simplified picture we are going to minimize the 
energy of the system upon making some approximations. 
We use Newton’s gravitational law (instead of general rel¬ 
ativity) and we assume a uniform density of DM fermions 
with Yukawa interactions. The energy has three contribu¬ 
tions, i.e. kinetic energy, gravitational potential energy and 
Yukawa potential energy 


E — Eq + Akin + Ey ■ (28) 


The self gravity contribution to the energy is 


Eg 



G4 
r 3 


7rpr^ • Airpr^dr 


3 Gm2 
5 R 


(29) 


where N is the total number of particles and R the radius of 
the star. The kinetic energy is found by multiplying Eq. (4) 
by the volume 


= (30) 

3 

where pf is the Fermi momentum of the particles. In the 
relativistic and non-relativistic limits, the kinetic energy 
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assumes simple polynomial forms in 1/i?: 


3. Relativistic, short range: pp ^ and /iR ^ 1: 


^^kin 


m^N + 


9s f Qtt 
Stt \2gs 


9s / Qtt \ 

IStt \2gJ 

^^4/3 ^4/3 


5/3 


N5/3 

m^R'^ 


In our approximation we define a radius 
PF(RrGi) = rn^, 


if pf ^ m^, 


if pf » m^. 

(31) 

i?rel such that 


_ f9n\ jVV3 

UsJ ’ 


(32) 


When R <C i?rei the system is relativistic, while in the 
opposite limit (i? i?rei) the system can be treated as 
non-relativistic. 

The Yukawa energy for the entire system of N particles is 
more complicated to derive than the other terms in the en¬ 
ergy. Starting from the potential between two particles in 
Eq. (8), one can find the potential by integrating the con¬ 
tributions of the shells of a homogeneous sphere of charges 
[81]. The final interaction energy reads 


Ey = -f 3 - 3(1 -f . 

^ 

This Yukawa potential energy also assumes simple polyno¬ 
mial forms in 1/i? in the long and short range limits where 
pi? ^ 1 and pi? 1 respectively 


^3aiV^ 
5 i? 

3 aN'^ 
^2p2i?3 


if pi? <C 1, 
if pi? ^ 1. 


(34) 


E{N, R) 


3 Gm2 iV2 / 97r \ ^ 3 aN^ 

5 R i? ^ 2 g?R^ ■ 

(37) 


4. Relativistic, long range: pp S> and pi? <C 1: 


E[N, R) 


3 Gm2 iV2 / Qtt \ iV4/3 3 aiV^ 

5 R i? ^ 5“^’ 

(38) 


We would like to estimate the condition for gravitational 
collapse of the objects, thus determining an upper mass 
limit for the asymmetric fermionic dark stars, similar to 
the Chandrasekhar limit in the case of white dwarfs. We 
assume that the star starts always from regime 1. However, 
depending on the parameters, the star might pass from 
different regimes before the collapse takes place. In practice 
we found two possibilities that take place most of the time. 
In the case of attractive interactions the star can move from 
regime 1 to 3 and collapse or from regime 1 to 3 to 4 and 
then collapse. For the repulsive potential the system goes 
directly from regime 1 to collapse. In the next subsections 
we will analyse each regime in greater detail. 


Regime 1: Non-relativistic, short range 


We rewrite Eq. (35) in terms of new constants H, B, C 
keeping the explicit dependence of N and R. 


E{N,R) = -A 


iV2 


+ B 


iV5/3 

i?2 


+ G 


/V2 

Ip' 


It is no surprise that the energy in the long range regime is 
independent of the mediator mass and proportional to 1 /R 
since it resembles the Coulomb potential. In the short range 
limit where pi? 3> 1, the exponential suppression of the 
potential is counterbalanced by the number of close neigh¬ 
bours, thus leading to a R~^ overall dependence. Given the 
above, we identify four distinct regimes that can be realised 
by the system: 

I. Non-relativistic, short range: pp ^ ^'^d pi? ^ 1: 

3 GmlN^ 

E(N,R) = -- ^ + m^N+ 

gs /97ry/^ iV5/3 3aN^ 

157 r\ 2 gJ m^R^ 2 p?R^' ^ ’ 


2. Non-relativistic, long range: pp and pi? <?/ I: 


E{N, i?) 


3 Gm2 iV2 
5 R 

gs / ^ 
IStt \2gs 


+ my^N+ 
5/3 ^5/3 

I 


3aN^ 
^ 5^~ 


(36) 


The derivative dE/dR = 0 gives two extrema 

For an attractive potential (G < 0) R^ is a minimum and 
i?_ a maximum. If one increases the value of N, there is a 
particular value where i?+ = i?_ and the system collapses 
because there is no stable solution. This value of N is 

/Mp^ 

V™x/ 
(40) 

When N = iV^ax, one can find that pi?+(Yniax) = 
3> I for the whole parameter space we 
examine here. Therefore when the system has accumulated 
Nmax particles, Yukawa forces are still short range. Since 
the energy has no minimum, i? keeps dropping until i?rei 
or the Schwarzschild radius i?s = 2,GNm^ is reached de¬ 
pending on which one is the larger. If i?rei > Rs, one can 
estimate pi?rBi (/Vmax) (2.6/ps)(I/x/a)Mpp2/m^ > 1. 
This means that for the parameter space considered here 
the system reaches i?rei before the Yukawa force becomes 


= - 


\ 

3^) 


ij'i 
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long range. Therefore the system will go to regime 3 or 4 
and not 2. If i?s > Rrd, the system collapses to a black 
hole before the particles become relativistic. In this case 
/ii?s ^ 1 for the parameter space we consider, and there¬ 
fore the system forms a black hole without reaching the 
long range limit of the Yukawa force. 

In the case where the Yukawa potential is repulsive 
(C > 0), R- < 0 (so it is unphysical) and R+ is a global 
minimum. The upper mass limit for this star can be found 
by setting i?+(iV) = Rs{N). There is always an appropri¬ 
ate N that satisfies this because R+{N) decreases with N 
while Rs{N) increases. For the parameter space we con¬ 
sider here, we found that the Schwarzschild radius is en¬ 
countered while still in the non-relativistic regime. The 
smallest possible value of i?+ (which occurs when TV —>■ oo) 
is Rs = \/‘iCjA. One can see that for the parameter space 
considered here, this line is crossed first by R^ and then by 
i?roi- It can be easily seen that this is also true even for a 
finite TV, thus the star collapses to a black hole before the 
constituents become relativistic. In addition, one can show 
that the aforementioned asymptotic value satishes [iR^ 3> 1 
and therefore the interactions are not long range. 


Regime 2: Non-relativistic, long range 


In the case of an attractive potential [C < 0), this ex¬ 
tremum is a maximum as long as TV is smaller than 


IVext — 



15 /MpY 

8 V 25s V / 


(44) 


For TV > TVext the energy is monotonically increasing as a 
function of R, and therefore the star collapses. For TV < 
TVext the energy increases between 0 and i?oxt and decreases 
from i?ext to oo. One can compare TV^ax from Eq. (40) and 
IVext 


IVmax _ 1 / /T \ ^ 

TVext " 125a3/2 Ux/ ■ 

If the star passes from regime 1 to regime 3, it collapses 
if IVmax > IVext- If On the other hand TV^ax < the 

collapse proceeds as long as i?rel(IVmax) < -Rext(IVmax)- 
In the case of repulsive Yukawa potential (C > 0), the 
potential is monotonically decreasing for TV < TVext- This 
means that the star remains in regime 1. Once TV > TVext 
there is a stable minimum at i?ext- However, as we pointed 
out in the discussion in ’’Regime 1”, the Schwarzschild ra¬ 
dius is met before this. 


We rewrite Eq. (36) with new constants A, B, C as 


Regime 4: Relativistic, long range 




In this case every term scales as \/R. Rewriting Eq. (38) 
in terms of new A, B, C constants we get 


The relevant minimum is found at 


R„ 


2B 

A-C 


TV-1/3 


Qtt \ 1 

2gs) 



(41) 

If the Yukawa potential is attractive {C < 0), there is al¬ 
ways a stable minimum. On the other hand in the case of 
repulsive potential C > 0, the strength of the coupling is 
bound by a < m^/TVfp, if a stable minimum is to exist at 
finite positive R. In practice for the considered parame¬ 
ter space, we found that no dark star passes through this 
regime. 


Tv2 TV^/3 tv^ 
E{N,R) = -A—+B—+C—. 


(46) 


■The critical number of particles is 


Nrrit — 


B 

A-c) 



Mv 


Ta 


-3/2 


(47) 

For TV > TVcrit the star collapses. In the case of repulsive 
interactions with a > m^/TVfp, no collapse can take place 
for any value of TV. As we mentioned in “Regime 1”, in case 
of repulsive interactions we never enter this regime because 
the particle never becomes relativistic. 


Regime 3: Relativistic, short range 

Once again we rewrite Eq. (37) in terms of new constants 
A, B, C as 


E{N, R) 



+ B 


TV4/3 

R 


TV2 


(42) 


VI. RESULTS 


In this section we present the full relativistic results af¬ 
ter solving numerically the system of Eqs. (23), (17), (11) 
and (12). The algorithm we use to solve the relativistic 
hydrostatic equilibrium is the following: 


and find the extremum at 


R 


ext — 




SC 


(43) 


1. Set initial conditions Mq = TW(r = 0) = 0 and 
Po = P{f = 0) = P{xo), with Xq = Xc being the 
relativity measure in the center of the star defined 
below Eq. (7). 


A-BN-‘^/3' 

















2. Integrate one step of Eq. (17) to get Mi = Mo+/ dM 
using the equation of state for pQ = p{xq). Then inte¬ 
grate one step of Eq. (23) to obtain Pi = Pq + J dP. 
For this value of Pi , one can obtain the corresponding 
xi from the equation of state (11). 

3. Repeat i times the above step to obtain Mi and Xi. 

4. Identify the R where P{R) = 0. This defines the 
radius of the star. Correspondingly the mass of the 
star is M{R). 

For each set of DM parameters we find the Chandrasekhar 
mass by scanning over x^ and identifying the largest total 
mass. 


VI.1. Mass-Radius relations 

We will present now the mass vs radius relations of 
the stable dark star configurations. In the left panels of 
Fig. 3 we show mass vs radius relations for the three 
generic cases: repulsive interactions (upper panel), no DM 
self-interactions (middle panel) and attractive interactions 
(bottom panel). We have chosen a coupling a = 10“^, 
mediator mass ^ = 10 MeV and three different cases of 
DM mass of 10, 100 and 1 TeV. For each M{R) profile, 
we mark the upper stable mass (Chandrasekhar mass) by a 
circle. Note that star configurations with radii larger than 
the one that corresponds to the Chandrasekhar mass are 
stable, while configurations with smaller radii are unstable. 
We refer the reader to the discussion of stability in subsec¬ 
tion IV. 1 and to the example of Fig. 2 which shows that the 
Chandrasekhar mass is indeed the last stable configuration 
as density increases. Note also that the dashed lines repre¬ 
sent dark star configurations where the non-relativistic ver¬ 
sion of the Tolman-Oppenheimer-Volkoff equation is used. 
As it can be seen, general relativity effects can be quite 
significant. 

For each mass vs radius curve that we show in the left 
panels of Fig. 3, apart from the Chandrasekhar mass, we 
mark in addition another point by a diamond. In the re¬ 
spective right panels, we show the density profile of the 
dark star that corresponds to the diamond. In the sce¬ 
narios of repulsive and no self-interactions, the equilibrium 
solutions feature a spiral structure in the unstable region of 
the M{R) curves. This feature is absent in the attractive 
scenario since the pressure becomes negative for some x, 
and this x is encountered before the spiral appears. One 
can notice (as it is expected) that heavier DM particles 
form lighter and more compact dark stars. 

In Fig. 4 we show again mass vs radius relations and 
corresponding density profiles for the diamond points (as 
in Fig. 3) for three values of a = 10“^, 10“^ and 10“'*, for 
fixed = 100 GeV and p, = 10 MeV for repulsive and at¬ 
tractive interactions. Note from Eqs. (11) and (12) that the 
Yukawa contribution to the density and the pressure of the 
star is unchanged under the scaling p ^ qp and a —)■ q^a. 
Therefore it is sufficient to fix either a and p and scan the 


phase space of the other one. One can notice that increas¬ 
ing the coupling of repulsive Yukawa interactions leads to 
larger stars, although DM particles in these stars are non- 
relativistic {xc <C 1). Despite that, the general relativity 
effects are large. 

VI.2. Chandrasekhar Mass 

In Fig. 5 we show the maximum mass of dark stars 
(Chandrasekhar limit) for the DM self-interaction param¬ 
eter space shown in Fig. 1. A first expected observation 
is that repulsive self-interactions lead to heavier dark stars 
compared to attractive ones, since repulsive interactions 
add to Fermi pressure and therefore more massive configu¬ 
ration can be supported. In fact larger couplings of a (for 
repulsive interactions) correspond to heavier stars. This 
can be seen also in Fig. 4. On the contrary, larger Yukawa 
couplings for attractive interactions lead to smaller Chan¬ 
drasekhar limits. One can see that the difference in the 
mass of dark stars made of DM particles with attractive 
and repulsive interactions is large. The mass of the attrac¬ 
tive dark stars for the whole parameter space lies below 
the limits that can be imposed by gravitational leasing. 
For repulsive stars, the mass can be a significant fraction 
of a solar mass, making these stars more visible from the 
leasing point of view. 

VII. DISCUSSION & CONCLUSIONS 

In this paper we studied the possibility that fermionic 
asymmetric DM can form stars, with self-interactions that 
can solve well standing problems of the CCDM paradigm. 
We solved the Tolman-Oppenheimer-Volkoff equation and 
we derived mass vs radius relations for these stars for both 
attractive and repulsive DM self-interactions for a range of 
DM and mediator masses. We also derived upper mass lim¬ 
its for these objects and we studied the hydrostatic stability 
of the star prohles we examined. 

One issue of fundamental importance that we plan to ad¬ 
dress in future work is related to the question of how these 
stars form in the first place. There are several possibilities 
where such a scenario can be visualized. One possibility 
is the creation of high DM density regions due to adia¬ 
batic contraction, caused by baryons [82] (see also discus¬ 
sion in [83]). Another possibility is the existence of a sub¬ 
dominant strongly interacting DM component that consists 
the seed for the gravitational collapse of parts of the DM 
halo via a gravothermal mechanism [84]. In this case DM 
self-interactions can lead to a collapse due to the fact that 
DM-DM collisions can send one of the two particles to a 
deeper gravitational potential while the second one leaves 
the area subtracting thus energy from the system that con¬ 
tinues until an instability is established. As it has been 
discussed in [55], the gravothermal mechanism leads to ei¬ 
ther gravitational collapse if the DM-DM forces are short 
range, or to formation of binary systems in case of long 
range interactions. In our case, the DM-DM interactions 
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R [km] r [km] 

(a) M{R) for repulsive interactions (b) p(r) for repulsive interactions 




r[km] 


(c) M (R) in the absence of interactions 


(d) p(r) in the absence of interactions 



FIG. 3: In the left panels we show dark star mass vs radius relations with DM mass = 10 GeV (Green), 100 GeV 
(blue), 1 TeV (purple). Upper, middle and bottom panels correspond to repulsive, no-interactions and attractive 
interactions respectively. We have fixed = 10 MeV and a = 10“^. Solid curves represent full relativistic solutions 
while dashed curves represent Newtonian gravity ones. The circles represent the Ghandrasekhar masses and the 
diamonds represent stars with their density profiles plotted as a function of the radius in the corresponding right 
panels. In the red regions R < Rg. In the attractive interaction scenario, the Newtonian solutions lie on top of the 
relativistic ones. 
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(a) M{R) for repulsive interactions 


(b) p{r) for repulsive interactions 



FIG. 4: In the left panels we show dark star mass vs radius relations for repulsive (upper panel) and attractive (right 
panel) for three different values of a = 10“^ (green), a = 10“^ (blue), and a = 10“"^ (purple). We have set = 100 
GeV, and /r = 10 MeV. For the star configurations of the diamond points, we show the corresponding density profiles 
on the respective right panels. 


are short range because they are mediated by the massive 
particle (j), (or cj)^) and therefore the instability leads to 
formation of DM asymmetric stars. Another possibility is 
the capture of a significant amount of DM particles by a 
supermassive star. After the collapse of the star and the su¬ 
pernova explosion, and since DM particles cannot be blown 
away in significant amounts by the supernova, a pure DM 
star or a mixed star with significant amount of baryons can 
be formed [85]. We should emphasize here that all these 
possibilities do not lead to a total collapse of the whole 
DM population but rather of a small fraction of it. In view 
of this, one should not worry for star profiles we present 
here in the case of repulsive interactions that are within the 
range of —10 Mq constrained by the MACHO [48] 

and EROS microlensing observations [49], simply because 
these stars do not comprise the whole DM density. 

Another issue we would like to leave out for future work 


is the discovery signatures for this type of dark stars. Apart 
from gravitational lensing that can in principle discover ob¬ 
jects like this based on the spacetime distortion that their 
presence can cause, other types of direct signals can exist. 
If DM communicates with the Standard Model via some 
portal, e.g. kinetic mixing between the photon and a dark 
photon, faint photon luminosity should be expected from 
these stars. Additionally, as we pointed out, for a large 
range of our parameters, the radius of these stars can be 
significantly smaller than that of a regular neutron star. 
This means that asymmetric dark stars can rotate faster 
than regular neutron stars. Pure dark stars made of DM 
particles that interact with the Standard Model particles 
through some portal, or mixed stars can appear as fast 
rotating pulsars. Rotational frequencies below millisecond 
are hard to be explained by a regular neutron star. An 
example of such a case is XTE J1739-285, which allegedly 























Mch [Mg] Mch [Mo] [Mq] 
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(a) a = 10 ^ attractive 



[GeV] 

(b) a = 10~^ repulsive 




(c) a = 10 ^ attractive (d) q = 10 ^ repulsive 



FIG. 5: Chandrasekhar mass for dark stars as a function of the DM mass for the parameter space of DM 
self-interactions shown in Fig. 1. 
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rotates with a frequency of 1122Hz [86]. Any odd look¬ 
ing neutron star is a potential candidate for an asymmetric 
dark star. 

We should mention that for a dark star that has collapsed 
to a black hole with a mass below ^ 10“ evaporation 

of the black hole via Hawking radiation takes place within 
the age of the universe. Although additional constraints 
exist [87-89], there is still a possibility for observing the 
spectrum of an evaporating black hole in the sky. 

Finally we should mention (as one can easily estimate) 
that in order for a dark star to collide with the earth 


during earth’s lifetime of to ^ 4.5 x 10® years, the dark 
star mass must he M < np^voR^to ~ where 

= 0.3GeV/cm^ is the local DM density, i?® is the radius 
of the earth, and Vq = 220km/sec the velocity dispersion 
of DM. One can see from Fig. 5 that ~TeV DM with 
attractive interactions and a coupling of a = 10“'^ can in 
fact give dark stars below that mass, making it possible for 
such a dark star to have collided with the earth in the past. 

This work is supported by the Danish National Re¬ 
search Foundation, Grant No. DNRF90. 
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